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ABSTRACT 


A  variety  of  two-equation  turbulence  models  -  including  several  versions  of  the  K  —  e 
model  as  well  as  the  K  -  u  model  -  are  analyzed  critically  for  near  wall  turbulent  flows  from 
a  theoretical  and  computational  standpoint.  It  is  shown  that  the  K  —  e  model  has  two  major 
problems  associated  with  it:  the  lack  of  natural  boundary  conditions  for  the  dissipation  rate 
and  the  appearance  of  higher-order  correlations  in  the  balance  of  terms  for  the  dissipation 
rate  at  the  wall.  In  so  far  as  the  former  problem  is  concerned,  either  physically  inconsistent 
boundary  conditions  have  been  used  or  the  boundary  conditions  for  the  dissipation  rate 
have  been  tied  to  higher-order  derivatives  of  the  turbulent  kinetic  energy  which  leads  to 
numerical  stiffness.  The  K  —  u>  model  can  alleviate  these  problems  since  the  asymptotic 
behavior  of  lo  is  known  in  more  detail  and  since  its  near  wall  balance  involves  only  exact 
viscous  terms.  However,  the  modeled  form  of  the  ui  equation  that  is  used  in  the  literature 
is  incomplete  -  an  exact  viscous  term  is  missing  which  causes  the  model  to  behave  in  an 
asymptotically  inconsistent  manner.  By  including  this  viscous  term  -  and  by  introducing 
new  wall  damping  functions  with  improved  asymptotic  behavior  -  a  new  K  —  r  model  (where 
r  =  l/u  is  turbulent  time  scale)  is  developed.  It  is  demonstrated  that  this  new  model  is 
computationally  robust  and  yields  improved  predictions  for  turbulent  boundary  layers. 


‘This  research  was  supported  by  the  National  Aeronautics  rind  Space  Administration  under  NASA  Con¬ 
tract  No.  NASl-18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 
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1.  INTRODUCTION 


An  increasing  number  of  practical  engineering  calculations  of  turbulent  flows  have  been 
based  on  two-equation  turbulence  models.  For  many  technologically  important  turbulent 
flows,  two-equation  models  represent  a  nice  compromise  between  zero  or  one  equation  mod¬ 
els  and  second-order  closures  (the  former  models  tend  to  require  too  many  ad  hoc  empiricisms 
whereas  second-order  closure  models  can  be  overly  expensive  for  design  calculations).  The 
K  —  c  model  [1]  is  the  most  popular  two-equation  turbulence  model  in  use  today.  When 
utilized  in  conjunction  with  wall  functions,  the  K  —  e  model  is  reasonably  well-behaved  and 
has  been  applied  to  the  solution  of  a  variety  of  engineering  problems  with  a  moderate  amount 
of  success.  However,  many  important  technological  applications  require  the  integration  of 
turbulence  models  directly  to  a  solid  boundary,  particularly  in  problems  where  wall  trans¬ 
port  properties  are  needed  or  where  there  is  flow  separation.  The  problem  of  developing 
low-Reynolds-number  near  wall  corrections  to  the  K  —  £  model  that  can  be  robustly  and 
accurately  integrated  to  a  solid  boundary  remains  unresolved  so  that  models  along  alterna¬ 
tive  lines  continue  to  be  proposed  (see  Patel,  Rodi,  and  Scheuerer  [2]  for  a  recent  review). 
Most  of  these  near  wall  K  —  e  models  involve  an  excessive  amount  of  ad  hoc  empiricisms 
and  are  numerically  stiff  in  turbulent  boundary  layer  flows.  This  motivated  some  researchers 
to  pursue  alternative  two-equation  models  based  on  a  modeled  transport  equation  for  the 
turbulent  time  scale.  The  most  notable  example  is  the  K  —uj  model  of  Wilcox  and  co-workers 
[3,  4]  where  modeled  transport  equations  for  the  turbulent  kinetic  energy  K  and  reciprocal 
turbulent  time  scale  uj  are  solved.  There  is  considerable  evidence  that  the  K  —  uj  model  is 
more  computationally  robust  than  the  K  —  £  model  for  the  integration  of  turbulent  flows 
to  a  solid  boundary.  However,  the  K  —  uj  model  yields  solutions  for  the  turbulent  kinetic 
energy  that  are  asymptotically  inconsistent  near  a  solid  boundary  [4].  Hence,  there  is  the 
need  to  re-examine  this  problem  from  a  basic  theoretical  and  computational  standpoint. 
This  establishes  the  motivation  for  the  present  paper. 

In  this  paper,  the  near-wall  asymptotics  of  two-equation  turbulence  models  will  be  ex¬ 
amined  from  a  basic  theoretical  standpoint.  It  will  be  shown  that  the  K  —  £  model  has 
two  major  problems  associated  with  it.  The  first  arises  from  the  lack  of  natural  boundary 
conditions  for  the  turbulent  dissipation  rate  which  has  caused  modelers  to  use  a  variety  of 
derived  boundary  conditions  that  are  either  asymptotically  inconsistent  (e.g.,  the  boundary 
condition  of  vanishing  normal  derivative  of  dissipation)  or  numerically  stiff  (e.g.,  the  bound¬ 
ary  condition  that  ties  the  dissipation  to  higher-order  derivatives  of  the  turbulent  kinetic 
energy).  The  second  problem  -  which  can  be  the  source  of  substantial  inaccuracies  and 
numerical  stiffness  -  is  tied  to  the  fact  that  the  balance  of  terms  at  the  wall  in  the  modeled 
dissipation  rate  transport  equation  depends  on  higher-order  correlations  whose  models  have 
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considerable  uncertainties. 

It  will  be  demonstrated  that  both  of  these  problems  can  be  largely  alleviated  by  solving 
a  modeled  transport  equation  for  the  turbulent  time  scale  t  —  K/e  since:  (a)  near  the  wall, 
r  y2/2i/  which  provides  the  needed  natural  boundary  conditions,  and  (b)  the  balance  of 
terms  at  the  wall  in  the  modeled  transport  equation  for  t  involves  only  exact  viscous  terms. 
It  will  be  argued  that  these  features  are  primarily  responsible  for  the  more  computationally 
robust  performance  of  the  K  —  uj  model  of  Wilcox  and  co-workers  [3,  4].  However,  the  K  —  u 
model  yields  results  for  the  turbulent  kinetic  energy  -  as  well  as  other  turbulence  quantities 
-  that  are  asymptotically  inconsistent  (e.g.,  near  the  wall,  K  ~  y3-23  instead  of  the  expected 
K  ~  y2  behavior).  It  will  be  proven  that  this  problem  arises  due  to  the  fact  that  an  exact 
viscous  cross-diffusion  term  is  missing  in  the  modeled  u;-transport  equation.  A  new  K  —  r 
model  is  obtained  by  including  this  exact  viscous  term  and  by  substituting  improved  wall 
damping  functions  which  are  obtained  by  an  asymptotic  analysis  using  the  results  of  direct 
numerical  simulations  of  turbulent  channel  flow  (Mansour,  Kim,  and  Moin  [5]).  The  new 
model  will  be  tested  for  the  flat  plate  turbulent  boundary  layer  and  comparisons  will  be 
made  with  the  predictions  of  other  models  (i.e.,  the  K  —  e  models  of  Launder  and  Sharma 
[6]  and  Chien  [7]  as  well  as  the  K  —  ui  model  [4])  in  order  to  assess  its  performance. 


2.  NEAR  WALL  ASYMPTOTIC  ANALYSIS 


For  simplicity,  we  will  restrict  our  analysis  to  incompressible  turbulent  flows  (however, 
the  crucial  conclusions  that  will  be  drawn  carry  over  to  compressible  flows).  The  mean 
velocity  u  and  mean  pressure  p  are  solutions  of  the  Reynolds  averaged  Navier-Stokes  and 
continuity  equations  given  by 


dui  _  dul 

~dt  +  U’lhi  = 


-£■ + + p 

oxi  oxj 


duj 

dxi 


=  0 


(1) 

(2) 
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where  Tjj  =  —u\u'j  is  the  Reynolds  stress  tensor,  v  is  the  kinematic  viscosity,  and  the  usual 
Einstein  summation  convention  applies  to  repeated  indices.  We  will  consider  the  commonly  1 

used  two-equation  models  based  on  an  eddy  viscosity  where 


TH 


+  vt 


(d*  &h\ 

\dxj  dxi) 


(3) 


and 


(4) 
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given  that  K  =  is  the  turbulent  kinetic  energy,  s  =  vdu'Jdxjdu'Jdxj  is  the  turbulent 

dissipation  rate,  and  C p  is  a  dimensionless  constant  at  high  turbulence  Reynolds  numbers. 
In  two-equation  models,  transport  equations  are  solved  for  any  two  linearly  independent 
variables  constructed  from  K  and  e.  In  the  K  —  e  model,  modeled  transport  equations  for  K 
and  e  are  solved;  in  the  K  ~ui  model,  modeled  transport  equations  for  K  and  the  reciprocal 
turbulent  time  scale  u>  =  c/K  are  solved;  and  in  the  K  —  r  model,  modeled  transport 
equations  for  K  and  the  turbulent  time  scale  r  =  K/e  are  solved.  The  exact  transport 
equations  for  K  and  e  are  as  follows  [8]: 


DK 

Dt 


=  T, 


dui 


-  e  -  V  +  uW2K 


where  D/ Dt  =  d/dt  +  u  •  V. 
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are  turbulent  transport  terms,  and 
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e  dxj  dxj  dxk  dx idxkdxk  dxkdxjdxj  kdx3dxjdxk 


=  2 1/2 


(10) 


dxjdxk  dxjdxk 

are,  respectively,  the  production  and  destruction  of  dissipation  terms. 

The  Taylor  series  expansions  for  the  components  of  u'  =  (u',v',w')  are  as  follows  near  a 
wall: 

u  =  aiy  +  a2y2  -I-  •  •  •  (11) 

v' =  b2y2  +  b3y3  +  ■  ■  ■  (12) 

w'  =  cry  +  c2y2  +  •  •  •  (13) 

where  a<  =  a,i(x,z,t),bi  =  bi(x,z,t)  and  c*  =  a(x,z,t)  given  that  the  coordinate  y  is  normal 
to  the  wall  (later,  wall  coordinates  will  be  used  where  y+  —  yuT/ v  given  that  ur  is  the  shear 
velocity).  Of  course,  the  no  slip  condition  at  the  wall  eliminates  ao,bo,  and  c0  whereas  the 
continuity  equation  eliminates  bi  (c.f.,  Hinze  [8]).  It  is  a  straightforward  matter  to  show  that 
near  a  wall 

K  =  0(y2),  e  =  0(  1),  r  =  0(y2)  (14) 
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rtTi  - 

—  =  0(1),  u'2  =  0(y2),  u'v'  =  0(y3) 

(15) 

^  =  0(y4 ),  ^  =  0(y2),  V  =  0(y 3) 

(16) 

V=0(y),  Ve  =  0(y),  *e  =  0(  1) 

(17) 

Vt  =  0(1),  V2A  =  0(1),  V2e  =  0(1) 

(18) 

where  V  =  Ttjdut/dxj  is  the  turbulence  production. 

An  asymptotic  analysis  of  the  K  —  e  model  will  be  conducted  first.  In  the  K  —  e  model, 
the  eddy  viscosity  near  a  wall  is  taken  to  be  of  the  form 

K2 


(19) 


The  asymptotic  analysis  presented  in  this  section  indicates  that  /M  =  0(1 /y)  near  the  wall 
since,  due  to  (15),  ut  must  be  of  0(y3)  in  this  region.  Of  course,  sufficiently  far  from  the 
wall  fp  assumes  a  value  of  1.  ( C ^  is  a  constant  which  is  typically  taken  to  be  0.09).  The 
turbulent  transport  term  T>  in  the  kinetic  energy  equation  (5)  is  modeled  using  a  gradient 
transport  hypothesis: 

f)  (20) 

where  aK  is  a  constant.  From  (14),  (17),  and  (19),  it  is  clear  that  this  model  is  not  asymp¬ 
totically  consistent.  However,  V  consists  of  two  parts  -  the  triple  velocity  term  and  the 
pressure  diffusion  term  -  as  given  by  (7).  Direct  numerical  simulations  of  the  Navier-Stokes 
equations  indicate  that 

d  ,1- 


a^(pu!)  <<: 

except  very  close  to  the  wall  (i.e.,  inside  of  y+  =  2;  see  Mansour,  Kim,  and  Moin  [5])  and  in 
this  region  V  is  negligible  in  comparison  to  the  dissipation  rate  and  the  viscous  diffusion  of 
the  turbulent  kinetic  energy.  Hence,  if  we  approximate  T>  by 


d  ,1- 


v  = 


(21) 


then  the  gradient  transport  model  (20)  is  asymptotically  consistent  since  the  right-hand-side 
of  (20)  and  (21)  are  both  of  0(y3)  as  the  wall  is  approached.  Hence,  it  would  appear  that 
the  asymptotic  errors  introduced  by  the  use  of  (20)  in  the  K  -  e  model  are  probably  not 
that  significant. 

The  turbulent  transport  term  Ve  in  the  dissipation  rate  transport  equation  is  also  modeled 
by  a  gradient  transport  hypothesis: 


vt  dc  ' 
ac  dxx 


(22) 


4 


(where  ae  is  a  constant)  in  the  K  —  e  model.  This  model  is  not  asymptotically  consistent  since 
T>s  —  0(1)  near  a  wall  while  the  right-hand-side  of  (21)  is  0{y2).  However,  this  inconsistency 
is  probably  not  of  great  consequence  since  both  T>t  and  $e  are  of  0(1)  near  a  wall  but  direct 
numerical  simulations  of  turbulent  channel  flow  indicate  that  Vc  <C  (f>e  (c.f.  Mansour,  Kim, 
and  Moin  [5]). 

The  production  of  dissipation  Ve  and  the  destruction  of  dissipation  $e  are  modeled  as 
follows 

p,  =  c>i/'7fT')^r  (23) 

=  C,ihj  (24) 

in  the  K  —  e  model  where  the  wall  damping  functions  fy ,  f2  — >  1  away  from  the  wall.  It 
is  clear  from  (14)  -  (17)  that  these  models  are  asymptotically  consistent  if  fy  =  0(1)  and 
f%  =  0(y2)  near  a  wall.  It  thus  follows  that  the  K  —  e  model  will  generate  solutions  for 
K,e  and  u'v'  that  are  asymptotically  consistent  if  the  damping  functions  /M  =  0(1  /y)  and 
f2  =  0(y2)  near  a  wall  with  fy  =  1. 

While  the  K  —  e  model  can  be  made  asymptotically  consistent  in  near  wall  turbulent 
flows  by  the  introduction  of  only  two  wall  damping  functions  -  namely,  f^  =  0(l/y)  and 
fi  =  0(y2)  -  there  are  still  some  other  major  problems  that  need  to  be  discussed.  There 
are  no  natural  boundary  conditions  on  e;  consequently,  boundary  conditions  must  be  either 
derived  or  postulated.  One  of  the  commonly  used  derived  boundary  conditions  is 


d2K 

V  »  =  £ 


dy2 


(25) 


at  the  wall  which  is  a  rigorous  consequence  of  the  exact  transport  equation  for  K .  Equation 
(25)  requires  information  at  the  wall  on  the  second-order  derivative  of  the  turbulent  kinetic 
energy  -  a  feature  that  can  lead  to  considerable  numerical  stiffness  [2].  Some  of  the  stiffness 
can  be  alleviated  by  utilizing  the  alternative  version  of  (25): 


-ft?)  - 


(26) 


at  the  wall.  However,  even  (26)  can  give  rise  to  considerable  numerical  problems.  The 
Neumann  boundary  condition 

fir 

(27) 


l  =  ° 


has  been  used  in  a  variety  of  applications  of  the  K  —  e  model  (c.f.  Lam  and  Bremhorst 
[9]).  Although  (27)  is  more  computationally  robust,  it  is  completely  ad  hoc  with  no  solid 
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theoretical  or  experimental  justification.  In  fact,  recent  direct  numerical  simulations  of  the 
Navier-Stokes  equations  for  turbulent  channel  flow  indicate  that  [5] 


1  de+ 
e+  dy+ 


-0.25 


at  the  wall;  under  such  circumstances  the  use  of  (27)  could  lead  to  substantial  errors. 

The  other  major  problem  with  the  dissipation  rate  transport  equation  lies  in  the  balance 
of  terms  at  the  wall.  At  a  solid  boundary,  (6)  reduces  to 


d2e  de 

‘+  ’  +  Tt 


(28) 


For  a  fully- developed  turbulent  boundary  layer,  de/dt  =  0  and  T>e  •<  3>e  as  discussed  earlier; 
hence,  (28)  simplifies  to 

(29) 


d2 £  . 

v - = 

dy 2  ‘ 


Both  (28)  and  (29)  have  a  major  deficiency:  the  balance  of  terms  at  the  wall  involves  higher- 
order  correlations.  This  puts  significant  pressure  on  the  accuracy  of  the  near  wall  modeling 
of  the  destruction  of  dissipation  term  that  can  further  exasperate  the  numerical  stiffness 
problem. 

On  the  other  hand,  the  turbulent  time  scale  r  =  KJe,  has  a  variety  of  natural  boundary 
conditions.  It  is  a  simple  matter  to  show  that  close  to  a  wall 


T  ~ 


r 

2  v 


and,  hence,  at  the  wall 


dr  d2r  1 

T  =  =  °>  -TT  =  -• 


dy 


dy2 


(30) 


(31) 


Equations  (30)  -  (31)  have  the  advantage  of  being  valid  for  any  near  wall  turbulence  where 
the  fluctuating  velocity  is  expandable  in  a  Taylor  series.  Furthermore,  the  balance  at  the 
wall  in  the  transport  equation  for  t  only  invokes  exact  viscous  terms.  This  can  be  seen 
directly  from  the  exact  r-transport  equation  which  takes  the  form 


Dt 

Dt 


T  dlii  T  _  T2  T2  T2 

~l~KV~  KVe  +  K*‘  +  KVe 


2v  dK  dr  2 v  dr  dr 


(32) 


K  dx{  dxi  t  dxx  dxi 

and,  hence,  at  the  wall  y  =  0,  the  leading  terms  are 


+  i/V2r 


2v  dK  dr  2v  dr  dr  d2r 
K  dy  dy  r  dy  dy  +  V  dy2 


(33) 
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Here,  each  term  on  the  left-hand-side  of  (33)  is  0(1).  The  balance  of  terms  in  (33)  is 
guaranteed  if  r  =  y2/ 2z/  near  the  wall.  It  therefore  appears  that  the  two  major  problems  with 
the  K  —  t  model  -  namely,  the  lack  of  natural  boundary  conditions  for  t  and  the  appearance 
of  higher-order  correlations  in  the  balance  of  terms  at  the  wall  -  can  be  overcome  by  the  use 
of  a  K  —  r  model. 

While  the  development  of  a  two-equation  turbulence  model  based  on  the  turbulent  time 
scale  has  been  discussed  in  the  literature  (c.f.  Reynolds  [10]  and  Bardina  [11]),  no  systematic 
study  of  such  models  has  been  conducted  for  near  wall  turbulence.  Only  the  K  —  u  model 
-  which  is  based  on  a  modeled  transport  equation  for  the  reciprocal  turbulent  time  scale 
<jJ  =  1/r  -  has  been  studied  in  these  flows  to  any  extent  (see  Wilcox  and  Traci  [3]  and 
Wilcox  [4]).  In  the  K  —  u>  model,  a  modeled  transport  equation  for  the  reciprocal  time  scale 
w  is  solved  which  is  of  the  form  [4] 


Duj  oj  dv,i  2  d  f  vt\  du 

~Dt  =  c“'KTiia^  “  °w2U  +  al,  [ v  +  Vj 


(34) 


where  ut  =  C^K/oj  and  Cui,CU2  and  are  constants  which  assume  the  values  of  5/9,  5/6, 
and  2,  respectively  (again,  C ^  =  0.09).  However,  we  will  now  show  that  (34)  is  inconsistent 
with  the  exact  transport  equation  for  w  near  a  wall:  an  exact  viscous  term  is  missing  and 
Cu 2  must  be  damped.  The  exact  transport  equation  for  u>  takes  the  form 


T r  Tf  Tr  ts  ^  ' 


u>T>  2u  du  dK 
+ 


+ 


(35) 


Dt  K  K  K  K  '  ~  '  K  '  K  dxidxi 
Hence,  it  is  clear  that  an  exact  viscous  cross-diffusion  term  -  given  by  (2 v / K)du / dxidK / dxi 
in  (35)  -  is  missing  in  the  modeled  w-transport  equation  of  Wilcox  and  co-workers.  From 
(35),  it  is  a  simple  matter  to  show  that  the  leading  terms  in  the  near  wall  balance  of  w  are 
as  follows: 

a2,., 

(36) 


2v  du>  dK  d2u>  2 

- 4-  V - f-  U  =  0 

K  dy  dy  dy2 


at  the  plane  solid  boundary  y  =  0.  Equation  (36)  is  consistent  with  asymptotic  solutions  for 
K  and  u  that  behave  correctly  near  the  wall,  i.e. , 

2v 


K  =  ay  ,  u;=-. 

y2 


(37) 


In  stark  contrast  to  (36),  the  K  —  u>  model  of  Wilcox  based  on  (34)  yields  the  balance  of 
terms 


d2u> 

dy2 


—  Cu  2<j)2  —  0 


(38) 


at  the  wall,  which  is  incompatible  with  (37).  Hence,  the  K  —  oj  model  of  Wilcox  yields 
asymptotically  inconsistent  solutions  in  near  wall  turbulence  (e.g.,  K  ~  y3,23;  see  Wilcox 

[4])- 
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The  K  —  uj  model  can  be  made  asymptotically  consistent  by  the  addition  of  the  viscous 
cross-diffusion  term 

2u  du  dK 
K  dxi  dxi 

and  by  decomposition  of  Cu 2  as  follows 

cu  2  =  C'u2-  1  (39) 

where  C'u2  must  be  damped  of  0(y2)  near  the  wall.  However,  we  feel  that  it  is  preferable 
to  derive  a  modeled  transport  equation  for  r  =  1/u;  since  r  is  not  singular  near  the  wall.  A 
new  K  —  t  model  will  be  derived  in  the  next  section  which  is  asymptotically  consistent. 


3.  A  NEW  K  —  t  MODEL 


The  exact  transport  equation  for  the  turbulent  time  scale  r  takes  the  form 

2i j  dr  dr 


Dt  t_  t  t2  t2  r2  2v  dK  dr 

Dt  K  K  K  K  K  K  dxi  dxi  r  ax{ 


+  i/V2r  (40) 


which  is  obtained  from  a  straightforward  combination  of  (5.)  and  (6).  Models  for  VC)  4>E,  T>t 
and  V  are  needed  for  closure.  The  production  of  dissipation  term  will  be  modeled  as  it  is  in 
the  K  —  £  model,  i.e., 

V  -  C  T  —  -  —t  — 

‘  ~  ,1K  ',9xj  ~  T  '’dx, 

where  c.  1  =  1.44.  As  mentioned  earlier,  this  model  is  asymptotically  consistent.  The 
destruction  of  dissipation  term,  will  also  be  modeled  similarly,  i.e., 

=  C.  (42) 

Here,  asymptotic  consistency  with  (17)  requires  that  f2  be  damped  of  0(y2)  near  a  wall.  We 
will  use  a  variation  of  the  form  for  Ce2  suggested  by  Hanjalic  and  Launder  [12]: 

Ce2  =  1.83[1  —  ^  exp(-Het/6)2]  (43) 

where  Ret  —  K2/ue  is  the  turbulence  Reynolds  number.  Here,  we  set  the  high  turbulence 
Reynolds  number  value  of  Ce2  =  1.83  since  it  yields  a  somewhat  more  accurate  value  for 
the  decay  rate  of  grid  turbulence  than  the  traditional  value  of  1.92  (c.f.  Reynolds  [10]).  An 
exponential  form  is  chosen  for  the  wall  damping  function  f2  as  follows 


fi  =  [1  -  exp(-y+/A2)]2 


(44) 
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(45) 


which  is  0(y2)  near  the  wall.  Since  at  the  wall 


d2e+ 
dy+ 2 


A2  can  be  evaluated  if  d2e+ /dy+2 ,  e+  and  K+  are  known.  By  using  the  wall  values  of  these 
quantities  from  direct  numerical  simulations  of  turbulent  channel  flow  [5]  we  obtain 


A2  =  4.9. 


The  resulting  model  for  <£e  is  quite  similar  to  that  proposed  recently  by  Myong  and  Kasagi 

[13]. 

The  turbulent  diffusion  term  for  r  is  defined  by 

Vr  =  K‘  ~  ic®'  (46) 

This  term  will  be  modeled  by  the  gradient  transport  hypothesis 

i/y  dr 
Or 2  dxx 

where  i/j  is  eddy  viscosity  and  an  and  a ^  are  turbulent  Prandtl  numbers.  In  deriving  (47) 
it  has  been  assumed  that  the  turbulent  transport  processes  parallel  the  molecular  ones  (i.e., 
each  turbulent  transport  term  is  coupled  with  a  molecular  diffusion  term  of  the  same  general 
form).  The  turbulent  Prandtl  numbers  for  the  last  two  terms  on  the  right-hand-side  of  (47) 
are  taken  to  be  equal  so  that  this  pairing  of  molecular  and  turbulent  diffusion  terms  is  also 
true  if  the  reciprocal  time  scale  uj  is  chosen  as  a  variable  instead  of  r  (the  choice  of  uj  as  a 
variable  instead  of  r  should  not  alter  the  basic  physics  of  the  model).  Of  course,  the  eddy 
viscosity  is  taken  to  be  of  the  form 


VT 


2 v?  dK  dr 
Kar\  dxi  dxi 


2z/y  dr  dr  d 

— i - 1 - 

tgt2  dxx  dxx  dxx 


VT  =  CJ^Kt  (48) 

where  C ^  =  0.09  and  is  a  wall  damping  function  which  is  0(  1/y)  near  the  wall.  By 
an  analysis  of  the  two  distinct  effects  of  low  turbulence  Reynolds  number  and  near  wall 
proximity,  Myong  and  Kasagi  [13]  proposed  the  model 

In  =  (1  +  3.45/y/Re^)[l  -  exp(— y+/70)].  (49) 

This  model  fits  the  experimental  data  [2]  reasonably  well  with  one  exception  -  it  is  asymptotes 
to  one  somewhat  too  slowly.  Hence,  we  will  consider  the  alternative  model 

U  =  (1  T  3.45/ \jHet)  tanh(y  +  /70)  (50) 
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since  the  hyperbolic  tangent  asymptotes  to  one  faster  by  the  necessary  amount. 

Now,  for  the  purposes  of  clarity,  we  will  summarize  the  K  —  r  model  derived  in  this 
section: 

2  (  dux  du,\ 

t„  =  --«u  +  „t(—  +  —  J  (51) 

vt  =  C^f^Kr  (52) 

D£_  dux  K  d\(  vr\dK] 

Dt  T'3dxj  t  +  dxi  [v  +  aK)  8xx\  ^ 

Dt  .  r,  \  r  ^Ui  tn  t  x\  2  /  Ur\dKdT 

Dt  -  {1  ~  ^K^dT,  +  {c‘2fl  ~ +  k  V  +  77J  dTiK 

(54) 

2  /  i/j  \  dr  dr  d  /  u j  \  dr 

r  \  ar2)  dxx  dxi  "**  dxx  V  ar2)  dxx 

In  (51)  -  (54),  C M  =  0.09,  =  1.44,  while  Ce2,  j M,  and  /2  are  given  by  (43),  (44)  and 

(50),  respectively.  However,  to  complete  the  model,  values  for  the  turbulent  Prandtl  numbers 
and  <7%  must  be  provided.  In  this  regard,  we  first  note  that  if 

&Ti  =  &T}  —  &K  —  &c 

then  the  modeled  r-transport  equation  (54)  is  equivalent  the  e  transport  equation 

De  e  dux  e2  d  \(  i/T\9el 

Dt  elKT'3dxj  Ce2^2K  +  dxi  [V  +  aj  dxx\  '  ^ 

Since  the  standard  e-transport  equation  (55)  is  known  to  perform  well  in  several  high- 
Reynolds-number  turbulent  flows,  it  is  reasonable  to  believe  that  an ,  ,  aK  and  ae  must 

assume  values  that  are  reasonably  close  to  one  another.  Furthermore,  for  local  equilibrium 
flows  with  zero  pressure  gradient  and  logarithmic  velocity  profile,  we  must  have  [2] 

c.  1  =  Ce2  ~  (K2/aeJc~„)  (56) 

where  k  =  0.4  is  the  Von  Karman  constant.  Hence,  for  the  values  of  Ce i,  Ce2  and  chosen 
herein,  it  follows  that 

cre  =  1.36 


and,  hence, 


<jr,  —  c-pj  —  <7/f  —  1.36. 
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It  should  be  noted  that  the  new  modeled  transport  equation  for  t  given  by  (54)  is 
equivalent  to  the  ^-transport  equation 


Du 

~Dt 


(C..  -  -  (C.,/2  -  l)w! 

/  vt  \  1  dK  du  d  (  vt  \  du 

V  crTl)  K  dx{  dxi  dxi  v  +  ari)  dx{ 


(58) 


which  differs  from  that  of  Wilcox  and  co-workers  by  the  presence  of  a  cross  diffusion  term 
and  by  the  damping  of  the  coefficient  of  u2  to  one  at  the  wall. 

Calculations  will  be  presented  in  the  next  section  using  the  common  value  of  1.36  for 
crn ,  a ^  and  aK  which  seems  to  be  adequate  for  the  present  study.  However,  future  research 
is  needed  to  optimize  these  constants  over  a  range  of  benchmark  turbulent  flows. 


4.  COMPARISON  OF  THE  MODELS 

Now,  the  performance  of  this  new  K  —  r  model  will  be  examined  for  the  flat  plate 
turbulent  boundary  layer  at  zero  pressure  gradient.  Comparisons  will  be  made  initially  with 
the  K  —  u  model  of  Wilcox  [4]  and  the  K  —  e  model  of  Launder  and  Sharma  [6]  (a  comparison 
with  the  K  —  e  model  of  Chien  [7]  will  be  made  later).  The  calculations  to  be  presented  were 
done  with  a  two-dimensional  boundary  layer  code  based  on  the  implicit  marching  scheme 
of  Edwards  et  al.  [14].  In  the  fully  developed  turbulent  regime,  approximately  100  grid 
points  were  used  in  the  direction  normal  to  the  wall  with  the  first  grid  point  at  y+  =  0.2. 
The  profiles  of  the  turbulent  fields  to  be  discussed  in  the  figures  are  for  a  Reynolds  number 
Ree  «  16,000  based  on  the  momentum  thickness  (this  will  allow  for  comparisons  with  the 
experimental  data  described  by  Patel,  Rodi,  and  Scheuerer  [2]  which  was  compiled  from  a 
variety  of  sources  including  Coles  [15]  and  Schubauer  [16]). 

In  Figure  1,  the  predictions  of  the  K  —  r,  K  —  u,  and  K  —  e  model  for  the  mean  velocity 
are  compared  with  experimental  data  [2].  It  is  clear  that  each  model  yields  a  logarithmic 
velocity  profile  for  30  <  y+  <  300  that  is  well  within  the  range  of  the  experimental  data. 
Furthermore,  each  model  correctly  yields  u+  =  y+  close  to  wall  (i.e.,  for  y+  <  5)  and 
predicts  the  deviations  from  the  law  of  the  wall  for  y+  >  1000.  In  Figure  2,  the  Reynolds 
shear  stress  predicted  by  these  three  models  is  shown.  The  predictions  of  the  various  models 
are  extremely  close  for  y+  >  10.  However,  for  y+  <  10  the  differences  between  the  model 
predictions  are  significant.  Among  these  models,  only  the  K  —  r  model  yields  a  profile  where 
u'v'  ~  y3  for  y+  <  10  as  indicated  by  experiments;  see  Figure  3  and  Patel  et  al.  [2].  In 
Figure  4,  the  predictions  of  the  K  —  r,  K  —  u>,  and  K  —  e  models  for  the  turbulent  kinetic 
energy  are  compared.  The  K  —  r  model  yields  a  peak  in  K+  of  approximately  4  which  is  well 
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within  the  range  of  the  experimental  data  [2]  and  the  results  of  direct  numerical  simulations 
for  turbulent  channel  flow  [5].  On  the  other  hand,  the  K  —  u>  model  -  as  well  as  the  K  —  e 
model  of  Launder  and  Sharma  -  appear  to  yield  peaks  in  the  turbulent  kinetic  energy  that 
are  rather  low.  The  turbulent  kinetic  energy  near  the  wall  is  shown  on  a  logarithmic  plot 
in  Figure  5.  Only  the  K  —  v  model  yields  K  ~  y2  for  the  entire  interval  0  <  y+  <  10;  it 
yields  the  proportionality  constant  ol+  =  0.05  -  a  result  that  is  well  within  the  range  of  the 
experimental  data. 

In  Figure  6,  the  profile  of  the  turbulent  dissipation  rate  predicted  by  the  K  —  r,  K  —  u> 
and  K  —  e  models  are  compared.  Although  the  results  are  fairly  close  for  t/+  >  20,  there 
are  some  significant  differences  close  to  the  wall.  The  K  —  t  model  yields  a  value  for  the 
turbulent  dissipation  rate  at  the  wall  of  £+  =0.1  which  is  quite  close  to  the  value  obtained 
from  experiments  [2].  Likewise,  the  peak  in  £+  is  quite  close  to  the  value  obtained  from 
experiments  [2].  In  contrast  to  these  results,  the  K  —  u  model  and  K  —  e  model  of  Launder 
and  Sharma  yield  values  for  the  wall  dissipation  that  are  substantially  too  small.  In 
Figure  7,  the  variation  of  with  y+  is  shown  for  these  three  models.  Only  the  K  —  r  model 
is  within  the  range  of  the  experimental  data  compiled  by  Patel,  Rodi,  and  Scheuerer  [2]. 

We  did  not  compare  directly  with  the  results  of  the  K  —  e  model  of  Chien  [7]  since  Chien’s 
model  was  calibrated  by  (and,  hen  »,  forced  into  agreement  with)  the  experimental  data  for 
the  flat  plate  turbulent  boundary  layer  at  zeTO  pressure  gradient.  However,  the  fact  that  this 
model  has  some  inconsistencies  can  be  seen  in  the  results  for  /M.  In  Figure  8,  a  comparison 
of  the  predictions  of  the  K  -  r  model  and  the  K  —  e  model  of  Chien  [7]  for  /M  is  shown.  It 
is  clear  from  this  figure  that  the  model  of  Chien  yields  values  for  /M  that  are  far  removed 
from  the  experimental  data.  Furthermore,  the  K  —  e  model  of  Chien  [7]  has  more  ad  hoc 
empiricisms  than  the  K  —  r  model  presented  herein. 

The  skin  friction  predicted  by  the  K  —  r  model  is  shown  as  a  function  of  the  coordinate 
x  along  the  plate  in  Figure  9.  It  is  clear  that  the  results  are  in  excellent  agreement  with  the 
experimental  data  [17].  In  Table  1,  the  fully-developed  skin  friction  and  wall  dissipation  rate 
are  tabulated  for  the  four  models  considered  in  this  study.  Only  the  K  —  r  model  and  K  —  e 
model  of  Chien  yield  results  that  are  within  the  range  of  the  experimental  data.  However,  it 
must  be  remembered  that  the  Chien  model  was  calibrated  by  forcing  it  into  agreement  with 
the  experimental  data  for  this  flow. 

The  more  desirable  features  of  r  as  a  variable  instead  of  e  can  be  seen  by  a  comparison 
of  Figure  10  with  Figure  6.  It  is  clear  from  Figure  10  that  the  turbulent  time  scale  varies 
much  more  smoothly  with  the  distance  from  the  wall;  its  first  derivative  with  respect  to  y 
does  not  change  sign. 
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5.  CONCLUSIONS 


A  basic  theoretical  and  computational  study  of  two-equation  models  for  near  wall  tur¬ 
bulent  flows  has  been  conducted.  The  major  findings  of  this  study  can  be  summarized  as 
follows: 

(1)  The  K  —  u>  model  of  Wilcox  and  co-workers  [3,  4]  is  missing  an  exact  viscous  cross 
diffusion  term.  Furthermore,  the  destruction  of  dissipation  term  is  not  properly  damped 
near  a  wall.  These  two  inconsistencies  give  rise  to  asymptotically  incorrect  solutions  for 
the  turbulent  kinetic  energy  ( K  ~  y3  23)  near  a  solid  boundary. 

(2)  The  K  —  e  model  can  be  made  asymptotically  consistent  by  the  satisfaction  of  two 
constraints:  the  coefficient  of  the  destruction  of  dissipation  term  must  be  damped  of 
0(y2)  near  a  wall,  and  the  coefficient  in  the  eddy  viscosity  must  be  damped  of  0(1 /y) 
near  a  wall.  Most  existing  corrections  to  the  K  —  e  model  yield  poor  results  in  near  wall 
turbulent  flows  due  to  the  violation  of  these  constraints. 

(3)  There  are  numerical  stiffness  problems  with  the  K  —  e  model  due  to  the  lack  of  natural 
boundary  conditions  for  the  dissipation  rate  and  the  fact  that  the  balance  of  terms  for 
the  dissipation  at  the  wall  involves  unknown  higher-order  correlations  which  need  to  be 
modeled.  These  problems  can,  to  a  large  extent,  be  overcome  by  the  use  of  the  turbulent 
time  scale  t  =  K/e  as  a  variable  since  r  has  natural  boundary  conditions  arising  from 
the  no-slip  condition  and  since  the  wall  balance  for  r  only  involves  exact  viscous  terms. 

(4)  A  new  K  —  t  model  was  developed  by  making  use  of  these  ideas  combined  with  a 
variation  of  wall  damping  functions  for  f2  and  that  were  recently  developed  by  Myong 
and  Kasagi  [13].  A  preliminary  test  of  this  K  —  r  model  for  the  turbulent  flat  plate 
boundary  layer  yielded  results  that  are  quite  encouraging.  However,  further  tests  and 
possible  refinements  are  required  before  more  definitive  conclusions  can  be  drawn. 

Future  research  will  be  directed  on  two  fronts.  The  K  —  r  model  will  be  subjected  to 
more  stringent  tests  involving  adverse  pressure  gradients  (with  possible  flow  separation)  and 
high  speed  compressible  flows.  The  turbulent  Prandtl  numbers  arl  and  ar2  will  be  optimized 
over  a  range  of  such  flows.  In  addition,  the  modeled  r-transport  equation  will  be  used  in 
conjunction  with  a  nonlinear  algebraic  stress  model  as  well  as  a  second-order  closure  model. 
One  of  the  major  goals  in  undertaking  this  research  was  to  ultimately  develop  second-order 
closure  models  that  can  be  integrated  directly  to  a  solid  boundary  in  complex  turbulent 
flows  that  involve  separation.  These  more  complex  issues  will  be  the  subject  of  our  ongoing 
research  effort  on  turbulence  modeling. 
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Model 

£\u 

Of 

K  —  T 

0.094 

0.00245 

K-e 

Launder-  Sharma 

0.021 

0.00238 

K-e 

Chien 

0.113 

0.00244 

K -u 

0 

0.00244 

Experiments  [2, 17] 

0.05-0.10 

0.00243 

Table  1.  Comparison  of  the  model  predictions  for  the  wall  dissipation  rate  ew  and  skin 
friction  Cj  ( x  =  4.987m,  Re  =  16,465,  and  Rex  =  1.156  x  107). 
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Figure  1.  Mean  velocity  profile  predictions  for  the  flat  plate  turbulent  boundary  layer  (Re$  « 

16, 000):  - K  -  t  model; - K  -  io  model;  •  •  •  K  -  t  model  of  Launder  and  Sharma 

[6];  O  experimental  data  [2]. 
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Figure  4.  Turbulent  kinetic  energy  profiles  (legend  the  same  as  in  Figure  1). 
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Figure  9.  Comparison  of  the  predictions  of  the  K  —  t  model  for  the  skin  friction  with 
experimental  data  [17]. 
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